The a — 7 transition of Cerium is entropy-driven 
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In this letter, we argue that entropic effects actually 
play the dominant role in the transition, at least in the 
temperature range 300K < T < T c . Using available ex- 
perimental data, we estimate the jump in entropy and 
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We emphasize, on the basis of experimental data and theoretical calculations, that the entropic 
stabilization of the 7-phase is the main driving force of the a-7 transition of cerium in a wide 
temperature range below the critical point. Using a formulation of the total energy as a functional 
of the local density and of the f-orbital local Green's functions, we perform dynamical mean-field 
theory calculations within a new implementation based on the multiple LMTO method, which allows 
to include semi-core states. Our results are consistent with the experimental energy differences and 
with the qualitative picture of an entropy-driven transition, while also confirming the appearance 
of a stabilization energy of the a phase as the quasiparticle Kondo resonance develops. 

PACS numbers: 71.27.+a, 71.30.+h, 71.15.Nc 

The a — 7 phase transition of cerium [1, 2] is a first- 
order isostructural transition, ending at a second-order 
critical point at T c ~ 600K. When temperature de- 
creases below T c , the volume change between the two 
phases increases [3, 4], reaching 15% at room temper- 
ature. The magnetic susceptibility follows Curie- Weiss 
behaviour in the (larger volume) 7- phase, and is Pauli- 
like in the (smaller-volume) a- phase. This is interpreted 
as 4f electrons being localized in the 7- phase, giving rise 
to local moments and contributing weakly to the elec- 
tronic bonding (hence the larger volume). In contrast, 
in the a- phase, the 4f electrons participate in both the 
bonding and the formation of quasiparticles. 

The detailed mechanism underlying the transition has 
been the subject of debate. In the Mott transition pic- 
ture [5], the focus is put on the 4f orbitals only, while 
the Kondo volume collapse (KVC) picture [6] empha- 
sizes the key role of the hybridisation between the 4f elec- 
trons and (spd-) conduction electron states which form 
broader bands. In this picture, the stronger hybridisa- 
tion of the low-volume a-phase leads to a high Kondo 
temperature and in turn to a screening of the 4f local 
moment, while the high- volume 7-phase has a low Kondo 
temperature, leading in practice to unscreened moments 
for T > Tj(. Photoemission experiments [7] demonstrate 
that both phases display Hubbard bands and hence are 
strongly correlated. In addition, a quasiparticle peak is 
seen in the a phase only. These observations are com- 
patible with both pictures. However, a recent theoretical 
calculation [8] of the optical spectrum, in connection with 
the experimental results of Ref. [9] , has emphasized the 
importance of hybridisation effects, in qualitative agree- 
ment with the KVC picture. In both the Mott and KVC 
pictures, the a- phase is stabilized by energetic effects 
(the f-electron kinetic energy in the Mott picture, or the 
Kondo screening energy in the KVC picture), while the 
7 phase is stabilized by its large spin-fluctuation entropy. 
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that TAS is always significantly larger than AE in this 
temperature range. The second purpose of this article is 
to examine whether this conclusion is consistent with to- 
tal energy calculations within the LDA+DMFT scheme, 
a combination of density-functional theory (DFT) within 
the local density approximation (LDA) with dynamical 
mean- field theory (DMFT). Recently, cerium has been 
the focus of pioneering theoretical work [8, 10-12] us- 
ing the LDA+DMFT approach. In Refs. [11], the total 
energy was studied and it was concluded that a nega- 
tive curvature effect is apparent already at elevated tem- 
peratures (T ~ 1600K), corresponding to the energetic 
stabilization of the cv-phase which was viewed as ulti- 
mately driving the transition. Here, we reconsider this 
issue within a new implementation of LDA+DMFT us- 
ing the multiple LMTO scheme, and basing our total- 
energy calculations on a functional of the local density 
and f-orbital Green's function. Using extensive Quan- 
tum Monte-Carlo calculations, we are able to reach tem- 
peratures lower than the experimental T c . Our results 
are consistent with the qualitative picture of an entropy- 
driven transition, and with the experimentally measured 
energy differences. 

The Clausius-Clapeyron relation dp/dT — AS/ AV re- 
lates the slope of the transition line to the jump of the 
entropy and unit-cell volume AV = V 1 — V a at the tran- 
sition. Furthermore, the continuity of the Gibbs free- 
energy yields the relation: AE — T AS + pAV = 0. Using 
available experimental data [3, 4] on dp/dT and AV, one 
can thus determine the three quantities AE, TAS and 
pAV, which are plotted on Fig. 1 as a function of temper- 
ature. As clear from this graph, the entropic term TAS 
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FIG. 1: Experimental variation of the entropy term TAS, en- 
ergy AE and pAV (obtained as described in the text from the 
data of [4] -circles- and [3] -squares-), between the 7 and the 
a phase across the transition line in the P-T phase diagram. 
Dotted lines are extrapolations based on exact limits (at T c 
all three terms vanish, and TAS = AE at p = 0). 



above room-temperature is of order 30 — 40meV, more 
than twice as large as the energy difference AE between 
the two phases (of order 10 — 20meV). The energetic sta- 
bilization of the a-phase results in AE — E 1 — E a > 0, 
but the difference in free-energy has the opposite sign: 
AF = F 1 - F a = AE - TAS = -pAV < precisely 
because the difference of entropy dominates over the en- 
ergy difference. We conclude from this analysis that en- 
tropic effects are essential to the physics of the a-j tran- 
sition, at least at room temperature and above. We note 
that the relative importance of spin and lattice contri- 
butions to the entropic stabilization of the 7-phase is 
currently under debate, with very different conclusions 
reached from experimental studies of pure cerium [13] 
and of the Ce .gTho.i alloy [14]. 

An accurate calculation of the electronic and lattice 
free-energy of cerium, from first principles, is a major 
challenge. While the calculation of the entropy is beyond 
the scope of this article, we focus in the following on the 
calculation of the total energy, at temperatures below the 
experimental T c , within the LDA+DMFT framework. 
This raises two important methodological questions. The 
first one is the proper choice of the valence states to be 
included in the starting LDA- hamiltonian. Indeed, it 
is mandatory to include semicore states (in particular 
5p states) in the valence when computing the energy, be- 
cause these states contribute significantly to the variation 
of the energy upon compression. On the other hand, we 
have found that, within an LMTO-ASA framework, it 
is crucial to include the 6p orbitals in the valence in or- 
der to obtain a proper band-structure. Since in standard 
implementations of the LMTO-ASA method, the simul- 
taneous inclusion of bp and 6p orbitals in the valence 
states is not possible, previous works using DMFT for 
Ce were either restricted to spectral properties only [10] 
or have treated the different terms in the expression of 



the total energy within different implementation of DFT 
[11]. Secondly, in view of the small energy differences 
between the two phases (on the scale of 10 to 20meV), 
a precise formulation of the total energy functional must 
be used. 

The starting Hamiltonian is constructed from an LDA 
calculation within the orthogonalized localized basis 
set of the multiple LMTO-ASA scheme [16], retain- 
ing 5s,5p,6s,6p,5d and 4/ states in the valence band. 
We neglect spin-orbit coupling which has little effect 
on LDA energies in Cerium. Many body terms act- 
ing on the f-orbitals are added to this hamiltonian, as 
well as a double-counting correction term (as in LDA+U 
schemes [17]), so that the many-body hamiltonian reads 
H = H KS + H V - H BC with: 



#KS= Z k LL>hfl(k)ci L C kL , 



(1) 



In this expression, /i^|,(k) denotes the Kohn-Sham 
(LDA) hamiltonian at a given k-point, expressed in an 
(orthogonalized) LMTO basis set Xlj with L = {Ima} 
and L' running over the full valence set. The Hubbard 
term is written in real space, with R denoting atomic po- 
sitions and a, b running only over the /-orbitals. We use 
the value of the Coulomb interaction U= 6 eV, computed 
by constrained LDA calculations in [2, 10]. 

In order to derive an expression for the total energy, 
we start from the ("spectral density") free-energy func- 
tional introduced by Kotliar and Savrasov [19], which de- 
pends on the total local electron density p(r) and the on- 
site Green's function in the correlated subset of orbitals: 
G RR (denoted below G ao for simplicity). The functional 
is constructed by introducing source terms, vks{ v ) ~ 
v c (y) (the difference of the Kohn Sham potential Vks 
and the crystal potential v c ), and A£ a (,(iw„), coupling 
to the density operators ip^ (r)ip(r) and to X^R.Xa( r ~ 

R)-0(r,r)?/' t (r',r')xb(i"' - R) = c q r(t)4 r (t'), respec- 
tively. The Luttinger-Ward [18] part of the functional 
is approximated by that of the on-site local many-body 
hamiltonian Hjj — i?DC- This yields: 

fi[p(r), G ab ; v K s(r), AS a b] lda+dmft 
= -itr In[«j„ + M + t^V 2 - v KS (r) - x*-A£ . X ] 

- J dr(«Ksr-« c )p(r)-tr[G.AE] 

+ I J dvdv'p{v)U{v-v')p{v') + E xc [p^)\ 

+ £(* imp [Gri-$Dc[G^ R ]) (2) 

R 

Minimization with respect to the sources gives a func- 
tional of the local Green function and the density only. 
Stationarity of this functional with respect to p(r) and 
G ab yields the basic equations of LDA+DMFT [20] , and 
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in particular, the self-consistency condition for the local 
Green's function: G a b(iu n ) = G(k, iu) n ) a t,. The full 
Green's function reads: Cr(k, iu n ) — (icu n +fj,)-l— h KS + 
V DC - t imp (iuj n ), with Eg^ p = 5§ imp /5G ab the local im- 
purity self-energy and V£q = 8$vc/3G a b- From (2), an 
expression of the total energy within LDA+DMFT can 
finally be obtained as: 



E = E.. 



DFT 



(H KS ) + (H v ) - E BC (3) 



data) for the a and the 7 phase, in comparison to ex- 
perimental spectra. As in previous LDA+DMFT stud- 
ies [10, 11], the quasiparticle peak is correctly described 
in the a phase, while Hubbard bands are present both 
in the a and the 7 phases. Their intensities are correct 
(although their positions are not very accurately repro- 
duced). These results give us confidence that the main 
physical features of both phases are correctly captured 
by our calculations. 



Note that, importantly, the total energy does not sim- 
ply reduce to the expectation value (H) of the many- 
body hamiltonian (1). In (3), -Edft is the expression 
of the energy within density-functional theory, 
is the sum of the occupied Kohn-Sham eigenvalues and 
(Hks) = tr [HksG]. Note that these last two terms do 
not cancel each other, since tr [HksG] is evaluated with 
the full Green's function including the self-energy, while 
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tr [HksGks]- Eq. (3) expresses that the lat- 
ter term has to be removed from Edft, in order to cor- 
rectly take into account the change of occupation of the 
Kohn Sham orbitals. The LDA+DMFT scheme should 
in principle be performed by imposing self-consistency 
not only on the DMFT quantities but also on the local 
density [15] (or equivalcntly on the LDA Hamiltonian), 
in such a way that p(r) = (r|G|r), i.c including the 
correlation-induced changes to the local density. How- 
ever, for simplicity and in order to compare to previous 
works [11], we present as a first step in this paper calcu- 
lations without full self-consistency on p(r) . The double- 
counting correction term is written in terms of the LDA 
occupancy of the f-orbital, as: Edc — UNf da (Ny da —l)/2, 
In order to solve the DMFT equations, we have used 
the Hirsch-Fye Quantum Monte Carlo algorithm, and 
studied the temperature range from 400K to 1600K. The 
number of sweeps was adjusted in order to obtain a pre- 
cision on the energy of order 20meV, a rather demanding 
requirement at the lowest temperature. The kinetic en- 
ergy (Hks) = En.k^stkJGfk,^) is computed in a 
direct manner, while the correlation energy (Hjj) is com- 
puted from the double occupancy. 
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FIG. 2: Experimental [7] and LDA+DMFT theoretical results 
for the PES and BIS spectra of a and 7 cerium 

First, we display on Fig. 2 the spectral functions (ob- 
tained by maximum-entropy continuation of our QMC 




FIG. 3: Symbols: internal energy (E) vs. volume curves for 
cerium, computed within LDA+DMFT for different temper- 
atures. Also shown is the free energy (F), calculated from the 
experimental pressure vs. volume curves (a: [13], b:[3]). The 
position of these experimental curves with respect to each 
other is arbitrary. Short vertical lines on these curves show 
the experimental volumes at each temperature. Arrows indi- 
cate the volume of each phase at room temperature, at the 
transition pressure. 

On Fig. 3, we display our results for the energy as 
a function of volume for three different temperatures 
(1600K, 800K and 400K). Statistical error bars of the 
QMC calculation are indicated on these plots. At 1600K 
we observe a smooth curve with a minimum located at 
31.0+0.5 A 3 . At 800K, the minimum is located at a 
somewhat lower volume and we note that the curva- 
ture decreases near the minimum of the curve and in 
particular for lower volumes. For 400K, this effect is 
strong enough to shift the minimum to 29+1 A 3 . Ac- 
cording to experimental results [3, 4], the volume of the 
a and 7 phases at 400K are 28.5+0.1 and 32.0+0.1 A 3 , 
and the difference of energy between these two phases 
is 13.5 ± 4meV (see Fig. 1). This value is quite consis- 
tent with our calculations, even though a precise theo- 
retical value would require to reduce the statistical error 
bars even more. Overall, we do not find evidence for 
a region of negative curvature in the energy versus vol- 
ume curve. However, because 400K is below the critical 
point, a double tangent should be present in the free- 
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energy vs. volume F(V). We have plotted in Fig. 3 
the experimental free-energy vs. volume curves deduced 
from recent pressure versus volume measurements at 413 
K [3] and 300K [3, 13], by integrating the equation of 
state: F(V) - F(V ) = - p(V')dV '. For volumes be- 
tween the equilibrium volumes of the a and 7 phases, 
(indicated by short vertical lines for each temperature), 
F(V) is taken to be the common tangent. Comparison 
of the theoretical energy to the experimental free-energy 
suggests that the entropic stabilization of the 7-phase is 
mainly responsible for the appearance of a region of neg- 
ative curvature in the free-energy F(V). Moreover the 
entropic stabilization TAS is of order 40meV at 400K, 
as seen from Figs. 1 and 3, much larger than AE. 

Although we do not find a double tangent in the en- 
ergy versus volume curve, we do observe a decrease of 
the curvature and the flattening of the volume depen- 
dance of the energy, as temperature is reduced. In or- 
der to understand its physical origin, we have plotted in 
Fig. 4, as a function of volume, the two contributions 
A = (Hks) — Sa£^ da (ie the correlation induced changes 
to the kinetic and hybridization energies) and B = (H\j) 
(the interaction energy among / orbitals). A negative 
curvature is clearly seen to develop in A as T is reduced, 
consistent with the observed development of the Kondo 
resonance (Fig. 2) and with the stabilization energy of 
the a phase as previously emphasized in reference [11]. 

Finally, the number of / electrons computed in DMFT 
is plotted in Fig. 5 as a function of volume. As expected, 
this number is very close to 1 in the localized 7 phase 
while it increases at lower volume due to hybridization ef- 
fects. In contrast to reference [11], we find a monotonous 
decrease of nt as volume is increased. One should keep 
in mind however that n/ depends on the basis set and 
the functional. 
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FIG. 4: Evolution of A = tr [H KS G] - E^ DA and B = (Hv) 
(note that, in this case, the scale is four time larger) as a 
function of volume, computed in DMFT for different temper- 
atures. 

In conclusion, we have revisited the problem of 
the volume-collapse transition of cerium, emphasizing 
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FIG. 5: Evolution of the number of f electrons in LDA and 
DMFT for different temperatures as a function of volume. 



that it is mainly entropy-driven. We have presented 
LDA+DMFT calculations of the total energy, obtained 
from a functional of the local density and local Green's 
function, within a new implementation based on the mul- 
tiple LMTO formalism. This allows us to include semi- 
core states and to calculate the Hamiltonian, the energy, 
as well as spectra within the same formalism. We confirm 
the development of a contribution to the kinetic and hy- 
bridisation energy stabilizing the a phase, as temperature 
is lowered and the Kondo quasiparticle resonance devel- 
ops, in qualitative agreement with the results of Ref. [11]. 
However, we find that the magnitude of this stabilization 
energy is too small to induce a pronounced negative cur- 
vature in the total energy curve, and that the transition 
is actually driven by entropy effects at least above room 
temperature. This is consistent with the experimental 
measurements of Drymiotis et al [21]. 
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orbit coupling. The conclusions of this recent work 
regarding cerium are qualitatively similar to those of 
Ref. [11]. 
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